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ABSTRACT 


We  describe  the  procedure  required  to  interface  the  FORTRAN  multiple  pre- 
cision package  of  Richard  P.  Brent  (as  described  in  ACM  Transactions  on  Mathe- 
matical Software,  March,  1978)  with  the  AUGMENT  precompiler  for  FORTRAN.  We 
also  indicate  the  method  of  using  the  multiple  precision  arithmetic  package  in 
conjunction  with  AUGMENT. 
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SIGNIFICANCE  AND  EXPLANATION 


In  some  applications,  it  is  necessary  to  use  higher  precision  than  is  af- 
forded by  standard  software.  The  multiple  precision  arithmetic  package  devel- 
oped by  Richard  P.  Brent  and  described  in  the  March,  1978  issue  of  ACM  Trans- 
actions on  Mathematical  Software  is  extremely  useful  in  such  cases. 

The  disadvantages  of  using  Brent's  package  directly  are  (1)  the  difficulty 
of  converting  existing  programs  to  make  use  of  the  multiple  precision  package, 
and  (2)  the  fact  that  in  order  to  write  a program  using  the  package,  one  must 
parse  the  arithmetic  expressions  oneself  and  write  the  program  as  a series  of 
calls  on  the  package  subroutines. 

The  AUGMENT  precompiler  for  FORTRAN,  developed  at  the  Mathematics  Research 
Center  by  F.  D.  Crary,  is  designed  to  simplify  the  use  of  packages  such  as 
Brent's.  In  this  report,  we  describe  the  necessary  Interface  to  enable  one  to 
use  Brent's  package  with  AUGMENT,  and  provide  instructions  for  its  use. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  sum- 
mary lies  with  MRC,  and  not  with  the  authors  of  this  report. 


AN  AUGMENT  INTERFACE  FOR  BRENT'S 
MULTIPLE  PRECISION  ARITHMETIC  PACKAGE 


Richard  P.  Brent,  Judith  A.  Hooper,  and  J.  M.  Yohe 


1.  Introduction: 

The  purpose  of  this  note  is  twofold:  first,  we  demonstrate  the  ease  with 
which  a well-designed  nonstandard  arithmetic  package  may  be  interfaced  with  the 
AUGMENT  precompiler  for  FORTRAN  [4];  second,  we  provide  an  interface  and  user 
instructions  to  enable  the  reader  to  use  Richard  P.  Brent's  FORTRAN  multiple 
precision  arithmetic  package  [1],  [2]  in  conjunction  with  AUGMENT.  This  makes 
the  use  of  Brent's  package  far  more  natural  and  convenient  than  its  use  without 
AUGMENT.  With  the  aid  of  AUGMENT,  the  user  declares  multiple  precision  varia- 
bles as  type  MULTIPLE,  and  then,  for  the  most  part,  simply  writes  the  program 
as  though  MULTIPLE  were  a standard  FORTRAN  data  type.  In  only  a few  instances 
must  the  user  write  explicit  calls  on  package  modules;  these  cases  will  be  dis- 
cussed later  in  the  paper. 


2.  Writing  the  interface: 

\ 

We  assune  that  the  reader  is  familiar  with  the  AUGMENT  precompiler,  at 
• least  to  the  extent  of  knowing  what  is  meant  by  such  terms  as  "supporting  pack- 

age" and  "description  deck".  This  degree  of  familiarity  may  be  gained  by  read- 
ing [4].  The  supporting  package  to  be  interfaced  with  AUGMENT  is  the  FORTRAN 
multiple  precision  arithmetic  package  described  by  Brent  in  [1]  and  [2].  This 
is  a collection  of  portable  subroutines  which  performs  not  only  basic  arithme- 
tic operations,  but  also  all  of  the  ANSI  standard  mathematical  functions  and 
many  nonstandard  ones,  in  multiple  precision.  The  precision  of  the  package  is 
governed  entirely  by  the  user  at  run  time,  and  may  even  be  changed  during  the 
course  of  a computation,  provided  the  dimensions  of  the  arrays  reserved  for  the 
multiple  precision  numbers  are  not  exceeded. 

■ 

In  interfacing  this  or  any  package  with  AUGMENT,  we  must  specify  the 
amoint  of  storage  to  be  allocated  to  each  variable.  This  will  place  an  upper 
limit  on  the  operating  precision  of  the  multiple  precision  arithmetic  package, 

I although  nothing  prevents  one  from  using  a lower  precision  in  computations. 

Increasing  the  precision  beyond  that  provided  in  this  standard  interface  is  not 
| [ difficult;  we  address  this  question  later. 

The  first  step  in  Interfacing  the  package  was  to  prepare  the  AUGMENT  de- 
scription decks.  The  multiple  precision  arithmetic  package  (although  not  de- 
i I signed  with  an  AUGMENT  interface  in  mind)  was  extremely  compatible  with  AUG- 
MENT: most  of  the  multiple  precision  routines  were  cast  as  subroutines,  with 


the  numbers  of  arguments  expected  by  AUGMENT  (and  in  the  order  expected);  near- 
ly all  of  the  manipulations  required  for  a complete  package  were  already  pro- 
vided (in  the  form  assumed  by  AUGMENT);  and  all  of  the  subroutines  in  the  pack- 
age bore  the  prefix  MP  in  their  names. 

The  preparation  of  the  description  deck  therefore  proceeded  easily;  we 
simply  went  down  the  list  of  routines  in  the  multiple  precision  arithmetic 
package,  associating  them  when  possible  with  standard  FORTRAN  operations  and 
functions.  When  such  a natural  association  was  not  possible,  we  assigned  func- 
tion names  (usually  obtained  by  dropping  the  prefix  'MP'  from  the  routine 
name) . The  description  of  each  routine  was  coded  as  per  the  instructions  in 
[5].  In  only  a few  cases  were  we  unable  to  do  this:  most  of  the  input/output 
routines  and  error  checking  routines  could  not  be  interfaced  with  AUGMENT  (they 
must  be  called  explicitly) , and  the  routines  which  provide  constants  needed 
special  attention,  as  we  shall  discuss  below.  Routines  which  did  not  conform 
to  the  usual  expectations  of  AUGMENT,  such  as  the  routine  to  add  the  quotient 
of  two  integers  to  a multiple  precision  number,  were  simply  described  as  func- 
tions. The  resulting  description  deck  is  shown  in  Appendix  B. 

The  routines  to  generate  constants  posed  a small  problem:  AUGMENT  assumes 
that  routines  will  have  at  least  one  argument  in  addition  to  the  result  (this 
might  be  regarded  as  a deficiency  in  AUGMENT),  and  these  routines  did  not.  We 
therefore  decided  to  write  a short  routine  to  interface  these  routines  with 
AUGMENT,  casting  it  as  a conversion  routine  which  "converts"  the  Hollerith  name 
of  the  desired  constant  to  the  value  of  the  constant.  This  routine  is  called 
with  the  Hollerith  name  of  the  constant  as  an  argument  (e.g.,  'PI'),  unpacks 
this  Hollerith  string,  determines  which  of  the  constant-generating  routines  to 
call,  and  returns  the  resulting  value.  Once  this  routine  was  written,  it 
seemed  logical  to  include  the  capability  of  run-time  conversion  of  numeric  con- 
stants, so  we  extended  the  routine  by  adding  a call  to  another  package  routine 
to  convert  the  (presumably  numeric)  Hollerith  string  to  multiple  precision  if 
it  did  not  match  the  name  of  any  of  the  "standard"  constants. 

We  also  wrote  six  trivial  logical  functions  to  allow  AUGMENT  to  deal  with 
the  six  logical  operators  in  the  context  of  multiple  precision  variables,  and 
some  other  routines  to  allow  the  user  to  inspect  and  modify  the  base,  number  of 
digits,  sign,  exponent,  and  digits  of  multiple  precision  numbers  without  need- 
ing to  know  the  details  of  the  implementation  of  the  package.  (These  should  be 
modified  only  with  extreme  care,  however.)  Finally,  we  added  some  input/output 
routines  which  are  simpler  to  use  with  the  AUGMENT  interface  than  those  origi- 
nally included  in  the  multiple  precision  arithmetic  package.  All  of  these  rou- 
tines were  extremely  straightforward  to  write  and  required  a total  of  about  120 
executable  statements.  A listing  of  them  is  given  in  Appendix  C. 

In  order  to  interface  the  PACK  and  UNPK  routines,  we  introduced  another 
data  type  called  MULTIPAK;  the  PACK  and  UNPK  routines  were  then  described  as 
conversions  between  types  MULTIPLE  and  MULTIPAK. 

The  entire  interface  was  written  in  less  than  a half-day;  the  most  time- 
consuming  task  was  revising  the  documentation  for  the  multiple  precision  pack- 
age! 
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3.  Use  of  the  package  via  AUGMENT: 


As  explained  in  [4],  the  use  of  a nonstandard  arithmetic  package  via  AUG- 
MENT is  extremely  simple.  The  majority  of  the  package  modules  are  invoked  au- 
tomatically by  AUGMENT,  the  exceptions  being  mainly  the  input/output  and  error 
handling  routines. 

To  use  the  package  through  AUGMENT,  the  user  declares  all  multiple  preci- 
sion variables  using  statements  of  the  form 

MULTIPLE  X,  Y(10),  Z 
or 

IMPLICIT  MULTIPLE  (A  - H,  0 - Z) 

(AUGMENT  accepts  type  declarations  via  IMPLICIT  statements,  whether  or  not  the 
FORTRAN  compiler  does;  this  is  convenient  when  converting  a program  to  multiple 
precision.)  The  majority  of  the  program  is  then  written  Just  as  though  MULTI- 
PLE were  a standard  FORTRAN  data  type. 

If  it  is  desired  to  store  multiple  precision  variables  in  packed  form,  one 
would  declare  a 10  by  100  array  of  packed  variables  in  the  following  way: 

MULTIPAK  A(10, 100) 

Since  the  package  normally  operates  only  on  unpacked  variables,  any  packed 
variables  must  normally  be  converted  to  unpacked  format  before  use.  This  may 
be  accomplished  by  either  of  two  methods: 

X s A(I,  J)  (normal  replacement  statement) 

CTM(A(I,  J))  (conversion  function) 

Packed  variables  should  not  normally  be  used  directly  in  arithmetic  ex- 
pressions, since  AUGMENT  will  not  generate  the  appropriate  conversion  in  all 
cases.  Packed  variables  may  be  used  in  certain  expressions;  for  example,  if  A 
and  C are  type  MULTIPLE  and  B is  type  MULTIPAK,  the  expression 

A = B • C 

will  work  properly.  However,  the  expression 
A = EXP(B) 

will  not  work;  it  must  be  written  as 
A s EXP(CTM(B) ) . 

The  user  may  elect  to  try  mixed  mode  expressions  of  other  kinds;  the  worst 
that  can  happen  is  that  the  linkage  editor  will  discover  that  AUGMENT  has  gen- 
erated a call  on  a nonexistent  routine. 

Constants  may  be  introduced  into  the  program  by  statements  of  the  follow- 
ing types: 

PI  * 'PI* 

X « ».1$» 
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The  dollar  sign  on  the  second  Hollerith  literal  Is  a sentinel  to  let  the 
Holler! th-unpacking  routine  know  when  it  has  reached  the  end  of  the  literal. 
If  the  compiler  generates  a sentinel,  and  if  the  unpacking  routine  recognizes 
it,  the  terminal  is  unnecessary.  (Note  that  the  Hollerith-unpacking  rou- 
tine is  NOT  portable;  it  will  need  to  be  rewritten  for  each  new  system.  The 
ones  shown  in  Appendix  C are  for  UNIVAC  1100  FORTRAN  V,  UNIVAC  1100  ASCII  FOR- 
TRAN, and  IBM  360  FORTRAN  G or  H,  respectively.) 

The  user  must  still  set  the  various  parameters  for  the  package,  as  ex- 
plained in  [1]  and  [33.  Care  must  be  exercised  to  ensure  that  the  dimensions 
of  the  multiple  precision  variables  communicated  to  the  package  are  no  greater 
than  those  used  by  AUOIENT  in  assigning  space  to  the  variables.  One  method  of 
setting  the  parameters  is  by  including  the  following  statements  in  the  main 
program,  before  any  (other)  executable  statements: 

COMMON  IDUMMY(K)  (where  K = MAXR  + 5) 

CALL  MPSET  (LUN,  NDIGIT,  N,  MAXR) 

where  LUN  is  the  logical  unit  nunber  for  output  (usually  6);  NDIGIT  is  the  num- 
ber of  decimal  digits  of  precision  desired;  N is  the  number  of  storage  loca- 
tions required  for  each  multiple  precision  variable  (this  must  not  exceed  the 
number  given  in  Line  23  of  the  description  deck  — 12  in  the  deck  shown  in  Ap- 
pendix B);  and  MAXR  is  the  length  of  the  working  space  array  as  described  in 
[1].  Of  course,  the  user  may  also  set  these  parameters  directly  as  described 
in  [33*  but  in  that  case,  care  must  be  exercised  not  to  exceed  the  number  of 
locations  assigned  to  variables  by  AUGMENT. 

Another  way  of  setting  these  parameters  to  default  values  is  to  include 
the  statement 

INITIALIZE  MP 

in  the  type  declarations.  This  causes  AUGMENT  to  generate  a call  on  the  rou- 
tine MPINIT,  which  then  sets  the  parameters  to  values  fixed  in  the  MPINIT  sub- 
routine. Of  course,  changes  in  the  dimensions  in  the  description  deck  must  be 
accompanied  by  appropriate  changes  in  the  parameters  in  MPINIT  if  this  method 
is  to  be  used.  This  is  a bit  of  a cludge,  but  it  works,  provided  the  default 
values  are  what  one  really  wants. 

A third  way  of  providing  these  parameters  to  the  package  would  be  even 
easier,  but  would  require  some  modification  of  the  package.  If  all  occurrences 
of  blank  COMMON  were  changed  to  labeled  COMMON  (e.g.,  00MM0N/MPC0M/) , the  pack- 
age parameters  could  then  be  set  via  a DATA  statement.  (This  was  not  done  in 
the  existing  package  because  of  a restriction  in  the  ANSI  (1966)  standards;  ac- 
cording to  these  standards,  labeled  COMMON  must  be  declared  consistently  in  all 
routines.)  The  setting  of  these  parameters  in  this  manner  would  obviate  the 
need  for  the  user  to  take  any  action  at  all;  however,  it  would  result  in  incom- 
patibility with  the  standard  (published)  version  of  the  multiple  precision 
arithmetic  package. 

The  maximum  precision  available  to  the  user  via  the  given  description  deck 
depends  on  the  word  length  of  the  host  computer;  on  the  UNIVAC  1110,  it  is  ap- 
proximately 43  digits.  The  value  of  MAXR  likewise  depends  on  the  characteris- 
tics of  the  host  machine  (and  on  the  modules  of  the  package  being  used  in  the 
program);  we  used  296  for  the  UNIVAC  1110  assuming  12  words  per  multiple  preci- 
sion variable. 


- 4 - 


i 


If  the  precision  provided  by  the  description  deck  is  not  sufficient  for 
the  user's  needs,  it  is  not  a difficult  task  to  increase  it;  one  merely  in- 
creases the  value  of  N given  in  Line  23  of  the  description  deck  to  accommodate 
the  desired  precision;  increases  the  value  in  Line  20  of  the  description  deck 
to  INT( (N  + 1)/2);  and  increases  MAXR  as  appropriate.  The  number  of  locations 
needed  for  the  work  space  array  will  depend  on  which  of  the  package  routines 
are  being  used;  the  amount  of  work  space  needed  for  each  routine  is  given  in 
[33.  The  most  space-consuming  routines  are  MPBESJ  and  MPLNGM . If  one  wishes 
to  avoid  the  pain  of  calculating  the  precise  requirements,  one  may  be  assured 
that  by  using 

MAXR  = max(T«*2  + 15*T  + 27,  14»T  + 156) 

where  T = N - 2,  enough  work  space  will  be  reserved  for  any  routine  in  the 
package.  These  considerations  are  discussed  in  greater  detail  in  [3]. 

Once  the  program  has  been  written,  the  following  runstream  will  invoke 
AUGMENT  and  cause  the  translated  program  to  be  written  on  logical  unit  20: 

(invoke  AUGMENT) 

(Description  Deck) 

•BEGIN 

(Source  Program) 

•END 

The  resulting  program  on  logical  unit  20  would  then  be  compiled  just  like 
any  other  FORTRAN  program;  the  compiled  program  would  then  be  linked  with  the 
multiple  precision  library  routines  and  executed. 

A complete  list  of  the  operations  and  functions  available  in  the  multiple 
precision  arithmetic  package,  together  with  the  manner  in  which  they  are  in- 
voked via  AUGMENT,  is  shown  in  Appendix  A. 


4.  Conclusion: 

We  have  demonstrated  the  method  of  interfacing  a supporting  package  with 
the  AUGMENT  precompiler  in  the  most  convincing  way  possible:  by  actually  doing 
it. 


The  interface  shown  in  this  paper  is  self-contained,  and  can  be  used  (with 
appropriate  modifications,  as  indicated  in  the  text)  with  Brent's  multiple  pre- 
cision arithmetic  package,  assuming  AUGMENT  is  available.  A revised  version  of 
the  multiple  precision  arithmetic  package,  incorporating  the  AUGMENT  interface 
routines,  is  available  from  the  first  author. 

Questions  may  be  addressed  to  the  authors. 
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APPENDIX  A 

OPERATIONS  IMPLEMENTED  IN  BRENT'S  MULTIPLE  PRECISION  PACKAGE 
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APPENDIX  A (Continued) 

OPERATIONS  IMPLEMENTED  IN  BRENT'S  MULTIPLE  PRECISION  PACK ACE 
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Check  legality  of  parameters  to  MP  package  - - MPCHKCIA,  IB) 

Handle  fatal  error  conditions  - - MPERR 

Handle  MP  overflow  - - MPOVFL(MR) 

Handle  MP  underflow  - - MPUNFL(MR) 


APPENDIX  A (Continued) 
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APPENDIX  B 


AUGMENT  DESCRIPTION  DECK  FOR  BRENT'S  MULTIPLE  PRECISION 
ARITHMETIC  PACKAGE 


•DESCRIBE  MULTIPAK 

COMMENT  AUGMENT  DESCRIPTION  DECK  FOR  THE  MULTIPLE-PRECISION 

ARITHMETIC  PACKAGE  OF  R.  P.  BRENT,  UNIVAC  1100  VERSION. 

THREE  TYPES  OF  VARIABLE  ARE  DEFINED  HERE  - 

MULTIPLE  (STANDARD  MULTIPLE-PRECISION  NUMBERS), 
MULTIPAK  (PACKED  MULTIPLE-PRECISION  NUMBERS),  AND 
INITIALIZE  (USED  ONLY  AS  A DEVICE  TO  PERSUADE 

AUGMENT  TO  INITIALIZE  THE  MP  PACKAGE). 
WORKING  SPACE  SHOULD  BE  ALLOCATED  AND  THE  MP  PACKAGE 
INITIALIZED  BY  THE  DECLARATION 
INITIALIZE  MP 
IN  THE  MAIN  PROGRAM. 

THIS  DESCRIPTION  DECK  ASSUMES  THAT  MULTIPLE  PRECISION  NUMBERS 
WILL  HAVE  NO  MORE  THAN  10  DIGITS  (BASE  65536)  FOR  A TOTAL 
PRECISION  NOT  EXCEEDING  ABOUT  43  DECIMAL  PLACES.  FOR  THIS, 
EACH  MP  NUMBER  REQUIRES  12  WORDS  (6  IN  PACKED  FORMAT). 

SEE  COMMENTS  IN  ROUTINE  MPINIT  FOR  THE  METHOD  OF  CHANGING 
THE  PRECISION  OR  ADAPTING  TO  A MACHINE  WITH  WORDLENGTH 
OTHER  THAN  36  BITS. 

DECLARE  INTEGERS),  KIND  SAFE  SUBROUTINE,  PREFIX  MPK 
SERVICE  COPY(STR) 

•DESCRIBE  MULTIPLE 

DECLARE  INTEGER(12),  KIND  SAFE  SUBROUTINE,  PREFIX  MP 
OPERATOR  + (.NULL  UNARY,  PRV,  $),  - (NEG,  UNARY), 

+ (ADD,  BINARY3,  PRV,  $,  $,  $,  COMM),  • (MUL) , 

- (SUB,,,,,,  NONCOMM),  / (DIV),  ••  (PWR2), 

♦ (ADDI,,,,  INTEGER),  • (MULI),  / (DIVI),  ••  (PWR), 

•EQ.  (EQ,  BINARY2,  PRV,  $,  LOGICAL,  COMM),  .NE.  (NE), 

.GE.  (GE,,,,,  NONCOMM),  .GT.  (GT),  .LE.  (LE),  .LT.  (LT) 

TEST  MPSIGA  (SIGA,  INTEGER) 

FIELD  SGN  (SIGA,  SIGB,  ($),  INTEGER), 

EXPON  (EXPA,  EXPB) , BASE  (BASA,  BASB) , NUMDIG  (DIGA,  DIGB) , 
MAXEXP  (MEXA,  MEXB) , DIGIT  (DGA,  DGB,  ($,  INTEGER)) 

FUNCTION  ABS  (ABS,  ($),  $),  ASIN  (ASIN),  ATAN  (ATAN),  CMF  ( CMF ) , 

CMIM  ( CMIM) , COS  (COS),  COSH  (COSH),  DAW  (DAW),  El  (El), 

ERF  (ERF),  ERFC  (ERFC) , EXP  (EXP),  EXP1  (EXP1 ) , FRAC  (CMF), 

GAM  (GAM),  INT  (CMIM),  LI  (LI),  LN  (LN),  LOG  (LN),  LNGM  (LNGM) 
LNGS  (LNGS),  LNS  (LNS) , REC  (REC),  SIN  (SIN),  SINH  (SINH), 
SQRT  (SQRT) , TAN  (TAN),  TANH  (TANH), 

ART1  (ART1,  (INTEGER)),  LN  (LNI),  LNI  (LNI),  LOG  (LNI), 

ZETA  (ZETA),  CAM  (CAM),  CAM  (CAM,  (HOLLERITH)), 

MAX  (MAX,  ($,  $)),  MIN  (MIN),  GCD  (GCDA), 

BESJ  (BESJ,  ($,  INTEGER)),  ROOT  (ROOT), 

MPINF  ( I NF( SUBROUTINE) , ( $ , INTEGER , INTEGER , HOLLERITH) , LOGICAL) , 
MPOUTF  (OUTF( SUBROUTINE)) , 

MPINF  (INF( SUBROUTINE),  ($,  INTEGER,  INTEGER,  INTEGER)), 

MPOUTF  (OUTF( SUBROUTINE)), 

COMP  (COMP,  ($,  $),  INTEGER),  CMPA  (CMPA), 

COMP  (CMPI,  ($,  INTEGER)),  COMP  (CMPR,  ($,  REAL)), 

ADDQ  ( ADDQ,  ($,  INTEGER,  INTEGER),  $) , MULQ  (MULQ) , 


MPA00010 
MPA00020 
MPA00030 
MPA00040 
MPA00050 
MPA00060 
MPA00070 
MPA00080 
MPA00090 
MPA00100 
MPA001 10 
MPA00120 
MPA00130 
MPA00140 
MPA00150 
MPA00160 
MPA00170 
MPA00180 
MPA00190 
MPA00200 
MPA00210 
MPA00220 
MPA00230 
MPA00240 
MPA00250 
MPA00260 
MPA00270 
MPA00280 
MPA00290 
MPA00300 
MPA00310 
MPA00320 
MPA00330 
MPA00340 
MPA00350 
MPA00360 
.MPA00370 
MPA00380 
MPA00390 
MPA00400 
MPA00410 
MPA00420 
MPA00430 
MPA00440 
MPA00450 
MPA00460 
MPA00470 
MPA00480 
MPA00490 
MPA00500 


QPWR  (QPWR,  (INTEGER,  INTEGER,  INTEGER,  INTEGER)), 

COM  (CQH,  (INTEGER,  INTEGER)),  CTM  (COM), 

GAM  (GAMQ) , GAMQ  (GAMQ) , 

BERN  (BERN,  (INTEGER,  INTEGER),  MULTIPAK) 

CONVERSION  CTM  (CDM,  DOUBLE  PRECISION,  $,  UPWARD), 

CTM  (CIM,  INTEGER),  CTM  (CRM,  REAL), 

CTM  (UNPK,  MULTIPAK),  CTM  (CAM,  HOLLERITH), 

CTD  (CMD( SUBROUTINE),  $,  DOUBLE  PRECISION,  DOWNWARD), 
CTI  ( CMI ( SUBROUT INE ) , , INTEGER), 

CTR  (CMR( SUBROUTINE),,  REAL),  CTP  (PACK,,  MULTIPAK) 
SERVICE  COPY  (STR) 

•DESCRIBE  INITIALIZE 

DECLARE  INTEGER ( 1 ) , KIND  SAFE  SUBROUTINE,  PREFIX  MPI 
SERVICE  COPY  (STR),  INITIAL  (NIT) 

COMMENT  END  OF  AUGMENT  DESCRIPTION  DECK  FOR  MP  PACKAGE 
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MPA00510 

MPA00520 

MPA00530 

MPA00540 

MPA00550 

MPA00560 

MPA00570 

MPA00580 

MPA00590 

MPA00600 

MPA00610 

MPA00620 

MPA00630 

MPA00640 

MPA00650 

MPA00660 


noooooooooooo 


APPENDIX  C 


AUGMENT  INTERFACE  ROUTINES  FOR  BRENT'S  MULTIPLE  PRECISION 

ARITHMETIC  PACKAGE 

C $$  «««  MPBASA  «*«»  MP009551 

FUNCTION  MPBASA  (X)  MP009553 

C RETURNS  THE  MP  BASE  (FIRST  WORD  IN  COMMON).  MP009555 

C X IS  A DUMMY  MP  ARGUMENT.  MP009557 

COMMON  B,  T,  M,  LUN,  MXR,  R MP009559 

INTEGER  B,  T,  R(1),  X(1)  MP009561 

MPBASA  = B MP009563 

RETURN  MP009565 

END  MP009567 

C $$  «»»*•  MPBASB  MP009571 

SUBROUTINE  MPBASB  (I,  X)  MP009573 

C SETS  THE  MP  BASE  (FIRST  WORD  OF  COMMON)  TO  I.  MP009575 

C I SHOULD  BE  AN  INTEGER  SUCH  THAT  I .GE.  2 MP009577 

C AND  (8*I*I-1 ) IS  REPRESENTABLE  AS  A SINGLE-PRECISION  INTEGER.  MP009579 

C X IS  A DUMMY  MP  ARGUMENT  (AUGMENT  EXPECTS  ONE).  MP009581 

COMMON  B,  T,  M,  LUN,  MXR,  R MP009583 

INTEGER  B,  T,  R(1),  X(1)  MP009585 

C SET  BASE  TO  I,  THEN  CHECK  VALIDITY  MP009587 

B = I MP009589 

CALL  MPCHK  (1,  4)  MP009591 

RETURN  MP009593 

END  MP009595 

C $$  •*••••  MPCAM  ••••••  MPO 12491 

SUBROUTINE  MPCAM  (A,  X)  MPO 1249 3 

CONVERTS  THE  HOLLERITH  STRING  A TO  AN  MP  NUMBER  X.  MP012495 

A CAN  BE  A STRING  OF  DIGITS  ACCEPTABLE  TO  ROUTINE  MPIN  MP012497 

AND  TERMINATED  BY  A DOLLAR  ($),  E.G.  7H-5.367A,  MP012499 

OR  ONE  OF  THE  FOLLOWING  SPECIAL  STRINGS  - MP012501 

EPS  (MP  MACHINE-PRECISION,  SEE  MPEPS),  MPO 12503 

EUL  (EULERS  CONSTANT  0.5772...,  SEE  MPEUL) , MP012505 

MAXR  (LARGEST  VALID  MP  NUMBER,  SEE  MPMAXR),  MP012507 

MINR  (SMALLEST  POSTIVE  MP  NUMBER,  SEE  MPMINR ) , MP012509 

PI  (PI  = 3.14...,  SEE  MPPI).  MP012511 

ONLY  THE  FIRST  TWO  CHARACTERS  OF  THESE  STRINGS  ARE  CHECKED.  MP012513 

SPACE  REQUIRED  IS  NO  MORE  THAN  5*T+L+14,  WHERE  L IS  THE  MP012515 

NUMBER  OF  CHARACTERS  IN  THE  STRING  A (EXCLUDING  $).  MP012517 

IF  SPACE  IS  LESS  3*T+L+11  THE  STRING  A WILL  EFFECTIVELY  BE  TRUNCATED  MP012519 
COMMON  B,  T,  M,  LUN,  MXR,  R MP012521 

INTEGER  B,  T,  R(1),  A(1),  X(1),  ERROR,  C(6),  D(2)  MP012523 

DATA  C(1)  /I HA/,  C(2)  /I HE/ , C(3)  /1HI/  MP012525 

DATA  C(4)  /1HM/,  C(5)  /1HP/,  C(6)  /1HU/  MP012527 

C UNPACK  FIRST  2 CHARACTERS  OF  A MP012529 

CALL  MPUPK  (A,  D,  2,  N)  MP012531 

IF  (N.NE.2)  GO  TO  10  MP012533 

C SET  X TO  ZERO  AFTER  SAVING  A(1)  IN  CASE  A AND  X COINCIDE  MP012535 

I * AO)  MP0 125 37 

* 0 MP0 12539 

C CHECK  FOR  SPECIAL  STRINGS  MPO 12541 
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IF  ((D(1).EQ.C(2)).AND.(D(2).EQ.C(5)))  CALL  MPEPS  (X)  MP012543 

IF  ((D(1).EQ.C(2)) . AND. (D(2) .EQ. C( 6) ) ) CALL  MPEUL  (X)  MP012545 

IF  ((D(1).EQ.C(4)) .AND. (D(2) .EQ. C( 1 ) ) ) CALL  MPMAXR  (X)  MP012547 

IF  ((D(1).EQ.C(4)) . AND. (D(2) .EQ.C(3) ) ) CALL  MPMINR  (X)  MP012549 

IF  ((D(1).EQ.C(5)).AND.(D(2).EQ.C(3)))  CALL  MPPI  (X)  MP012551 

C RETURN  IF  X NONZERO  (SO  ONE  OF  ABOVE  TESTS  SUCCEEDED)  MPO 12553 

IF  (X(I).NE.O)  RETURN  MP012555 

C RESTORE  A(1)  AND  UNPACK,  THEN  CALL  MPIN  TO  DECODE.  MPO 12557 

A(1)  = I MPO 12559 

10  12  s 3«T  + 12  MPO 12561 

CALL  MPUPK  (A,  R(I2),  MXR+1-I2,  N)  MP012563 

CALL  MPIN  (R(I2),  X,  N,  ERROR)  MP012565 

IF  ( ERROR. EQ.O)  RETURN  MPO 12567 

WRITE  (LUN,  20)  MPO 12569 

20  FORMAT  (53H  •**  ERROR  IN  HOLLERITH  CONSTANT  IN  CALL  TO  MPCAM  »••)  MPO 12571 

CALL  MPERR  MPO 12573 

RETURN  MPO 12575 

END  MPO 12577 


C $$  MPDGA 

FUNCTION  MPDGA  (X,  N) 

C RETURNS  THE  N-TH  DIGIT  OF  THE  MP  NUMBER  X FOR  1 .LE.  N .LE.  T. 
C RETURNS  ZERO  IF  X IS  ZERO  OR  N .LE.  0 OR  N .GT.  T. 

COMMON  B,  T,  M,  LUN,  MXR,  R 
INTEGER  B,  T,  R(1),  X(1) 

MPDGA  = 0 

IF  ((X(I).NE.O) .AND. (N.GT.O) .AND. (N. LE.T) ) MPDGA  = X(N+2) 

RETURN 

END 


MPO 19741 
MPO 19743 
MPO 19745 
MPO 19747 
MPO 19749 
MPO 19751 
MPO 19753 
MPO 19755 
MPO 19757 
MPO 1975 9 


C $$  MPDGB  •••••• 

SUBROUTINE  MPDGB  (I,  X,  N) 

C SETS  THE  N-TH  DIGIT  OF  THE  MP  NUMBER  X TO  I. 

C N MUST  BE  IN  THE  RANGE  1 .LE.  N .LE  T, 

C I MUST  BE  IN  THE  RANGE  0 .LE.  I .LT.  B 
C (AND  I .NE.  0 IF  N .EQ.  1). 

C THE  SIGN  AND  EXPONENT  OF  X ARE  UNCHANGED. 

COMMON  B,  T,  M,  LUN,  MXR,  R 
INTEGER  B,  T,  R(1),  X(1) 

IF  ((N. GT.O). AND. (N. LE.T))  GO  TO  20 
WRITE  (LUN,  10) 

10  FORMAT  (48H  •••  DIGIT  POSITION  ILLEGAL  IN  CALL  TO  MPDGB  •••) 
GO  TO  40 

20  IF  ( (I .GE.O) .AND. (I.LT.B) .AND. ( (I+N) .GT. 1 ) ) GO  TO  50 
WRITE  (LUN,  30) 

30  FORMAT  (45H  •••  DIGIT  VALUE  ILLEGAL  IN  CALL  TO  MPDGB  •••) 

40  CALL  MPERR 
RETURN 

50  X ( N+2 ) = I 
RETURN 
END 


C $$  «••••  MPDIGA  •••*•• 

FUNCTION  MPDIGA  (X) 

C RETURNS  THE  NUMBER  OF  MP  DIGITS  (SECOND  WORD  IN  COMMON). 
C X IS  A DUMMY  MP  ARGUMENT. 


MP019781 
MPO 19783 
MPO 19785 
MPO 19787 
MPO 19789 
MPO 19791 
MPO 1979 3 
MPO 19795 
MPO 19797 
MPO 19799 
MPO 19801 
MPO 1980 3 
MPO 19805 
MPO 19807 
MPO 19809 
MP019811 
MP019813 
MPO 198 15 
MP019817 
MP019819 
MPO 19821 

MPO 19841 
MPO 1984 3 
MPO 19845 
MPO 19847 


- 14  - 


W 


COMMON  B,  T,  M,  LUN,  MXB,  R MP019849 

INTEGER  B,  T,  R(1),  X(1)  MP019851 

MPDIGA  = T MPO 19853 

RETURN  MPO 19855 

END  MPO 19857 

C $$  •*•••«  MPDIGB  »•••••  MPO 19861 

SUBROUTINE  MPDIGB  (I,  X)  MP019863 

C SETS  THE  NUMBER  OF  MP  DIGITS  (SECOND  WORD  OF  COMMON)  TO  I.  MP019865 

C I SHOULD  BE  AN  INTEGER  SUCH  THAT  I .GE.  2 MPO 19867 

C X IS  A DUMMY  MP  ARGUMENT  (AUGMENT  EXPECTS  ONE).  MP019869 

C WARNING  •••  MP  NUMBERS  MUST  BE  DECLARED  AS  INTEGER  ARRAYS  OF  MPO 19871 

C •••  DIMENSION  AT  LEAST  1+2.  MPDIGB  DOES  NOT  CHECK  THIS.  MPO 1987 3 

COMMON  B,  T,  M,  LUN,  MXR,  R MP019875 

INTEGER  B,  T,  R(1),  X(1)  MP019877 

C SET  DIGITS  TO  I,  THEN  CHECK  VALIDITY  MP019879 

T = I MPO 19881 

CALL  MPCHK  (1,4)  MPOI9883 

RETURN  MPO 19885 

END  MP019887 

C $$  ••••••  MPEQ  »«•«  MP023221 

LOGICAL  FUNCTION  MPEQ  (X,  Y)  MP023223 

C RETURNS  LOGICAL  VALUE  OF  (X  .EQ.  Y)  FOR  MP  X AND  Y.  MP023225 

INTEGER  X(1),  Y ( 1 ) MP023227 

MPEQ  = (MPCOMP(X,Y)  .EQ.  0)  MP023229 

RETURN  MP023231 

END  MP023233 

C $$  ••••••  MPEXPA  MP027271 

FUNCTION  MPEXPA  (X)  MP027273 

C RETURNS  THE  EXPONENT  OF  THE  MP  NUMBER  X MP027275 

C (OR  LARGEST  NEGATIVE  EXPONENT  IF  X IS  ZERO).  MP027277 

COMMON  B,  T,  M,  LUN,  MXR,  R MP027279 

INTEGER  B,  T,  R(1),  X(2)  MP027281 

MPEXPA  = -M  MP027283 

C RETURN  -M  IF  X ZERO,  X(2)  OTHERWISE  MP027285 

IF  (X(I).NE.O)  MPEXPA  = X(2)  MP027287 

RETURN  MP027289 

END  MP027291 

C $$  ••••«  MPEXPB  ••••••  MP027311 

SUBROUTINE  MPEXPB  (I,  X)  MP027313 

C SETS  EXPONENT  OF  MP  NUMBER  X TO  I UNLESS  X IS  ZERO  MP027315 

C (WHEN  EXPONENT  IS  UNCHANGED).  MP027317 

C X MUST  BE  A VALID  MP  NUMBER  (EITHER  ZERO  OR  NORMALIZED).  MP027319 

COMMON  B,  T,  M,  LUN,  MXR,  R MP027321 

INTEGER  B,  T,  R(1),  X(3>  MP027323 

C RETURN  IF  X IS  ZERO  MP027325 

IF  (X(I).EQ.O)  RETURN  MP027327 

C CHECK  FOR  VALID  MP  SIGN  AND  LEADING  DIGIT  MP027329 

IF  ( ( IABS(X( 1 ) ) .LE. 1 ) .AND. (X(3) >GT.O) .AND. (X(3) . LT. B) ) MP027331 

$ GO  TO  20  MP027333 

WRITE  (LUN,  10)  MP027335 

10  FORMAT  (48H  •••  X NOT  VALID  MP  NUMBER  IN  CALL  TO  MPEXPB  •••)  MP027337 

CALL  MPERR  MP027339 
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X(1)  s o 
RETURN 

C SET  EXPONENT  OF  X TO  I 
20  X(2)  = I 

C CHECK  FOR  OVERFLOW  AND  UNDERFLOW 
IF  (I.GT.M)  CALL  MPOVFL  (X) 

IF  (I.LT.(-M))  CALL  MPUNFL  (X) 

RETURN 

END 


C $* 


«»•*•  MPGE 

LOGICAL  FUNCTION  MPGE  (X,  Y) 

C RETURNS  LOGICAL  VALUE  OF  (X  .GE.  Y)  FOR  MP  X AND  Y. 
INTEGER  X(1),  Y( 1 ) 

MPGE  = (MPCOMP(X.Y)  .GE.  0) 

RETURN 

END 


C $* 


**••••  MPGT  »»*»»* 

LOGICAL  FUNCTION  MPGT  (X,  Y) 

C RETURNS  LOGICAL  VALUE  OF  (X  .GT.  Y)  FOR  MP  X AND  Y. 
INTEGER  X(1),  Y ( 1 ) 

MPGT  = (MPCOMP(X.Y)  .GT.  0) 

RETURN 

END 


C $$ 


»••***  MPINF 

SUBROUTINE  MPINF  (X,  N,  UNIT,  IFORM,  ERR) 

READS  N Words  from  logical  UNIT  IABS(UNIT)  using  FORMAT  in  IFORM, 
THEN  CONVERTS  TO  MP  NUMBER  X USING  ROUTINE  MPIN. 

IFORM  SHOULD  CONTAIN  A FORMAT  WHICH  ALLOWS  FOR  READING  N WORDS 
IN  A1  FORMAT,  E.G.  6H(80A1) 

ERR  RETURNED  AS  TRUE  IF  MPIN  COULD  NOT  INTERPRET  INPUT  AS 
AN  MP  NUMBER  OR  IF  N NOT  POSITIVE,  OTHERWISE  FALSE. 

IF  ERR  IS  TRUE  THEN  X IS  RETURNED  AS  ZERO. 

SPACE  REQUIRED  3T+N+11. 

COMMON  B,  T,  M,  LUN,  MXR,  R 

INTEGER  B,  T,  R(1),  X(1),  UNIT,  IFORM(I) 

LOGICAL  ERR 

CHECK  THAT  ENOUGH  SPACE  AVAILABLE 
CALL  MPCHK  (3,  N+11) 

12  = 3*T  ♦ 12 

READ  N WORDS  UNDER  FORMAT  IFORM. 

CALL  MPIO  ( R ( 12) , N,  (-IABS(UNIT)) , IFORM,  ERR) 

X(1)  = 0 
RETURN  IF  ERROR 

IF  (ERR)  RETURN 
ELSE  CONVERT  TO  MP  NUMBER. 

CALL  MPIN  (R(I2),  X,  N,  IER) 

RETURN  ERROR  FUG  IF  MPIN  OBJECTED 
ERR  : (IER.NE.O) 

RETURN 

END 


C $$ 


MPINIT 


MP027341 

MP027343 

MP027345 

MP027347 

MP027349 

MP027351 

MP027353 

MP027355 

MP027357 

MP030521 
MP030523 
MPO 30525 
MP030527 
MP030529 
MP030531 
MP030533 

MP030541 

MP030543 

MP030545 

MP030547 

MP030549 

MP030551 

MP030553 

MP032761 
MP032763 
MP032765 
MP032767 
MP032769 
MP032771 
MP032773 
MP032775 
MP032777 
MP032779 
MP032781 
MP032783 
MP032785 
MP032787 
MP032789 
MP032791 
MP032793 
MP032795 
MP032797 
MP032799 
MP032801 
MP032803 
MP032805 
MP032807 
MP032809 
MPO 3281 1 
MP032813 

MP032821 
MPO 3282 3 


1 


SUBROUTINE  MPINIT  (X) 
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C DECLARES  BLANK  COMMON  (USED  BY  MP  PACKAGE)  AND 
C CALLS  MPSET  TO  INITIALIZE  PARAMETERS 
C THE  AUGMENT  DECLARATION 
C INITIALIZE  MP 

C CAUSES  A CALL  TO  MPINIT  TO  BE  GENERATED. 

C «*  ASSUMES  OUTPUT  UNIT  6,  43  DECIMAL  PLACES, 

C •••  10  MP  DIGITS,  SPACE  296  WORDS.  IF  THE  AUGMENT 
C •••  DESCRIPTION  DECK  IS  CHANGED  THIS  ROUTINE  SHOULD 
C •••  BE  CHANGED  ACCORDINGLY. 

COMMON  B,  T,  M,  LUN,  MXR,  R 
INTEGER  B,  T,  X( 1 ) 

C THE  STATEMENTS 

INTEGER  R(296) 

CALL  MPSET  (6,  43,  12,  296) 

C ARE  A SPECIAL  CASE  OF 
C INTEGER  R(MXR) 

C CALL  MPSET  (LUN,  IDECPL,  T+2,  MXR) 

C WHERE  LUN  IS  THE  LOGICAL  UNIT  FOR  OUTPUT, 

C IDECPL  IS  THE  EQUIVALENT  NUMBER  OF  DECIMAL  PLACES  REQUIRED, 

C T IS  THE  NUMBER  OF  MP  DIGITS,  AND 
C MXR  IS  THE  SIZE  OF  THE  WORKING  AREA  USED  BY  MP 
C (MXR  = MAX  (T«T+15*T+27,  14»T+156)  IS  SUFFICIENT). 

C TO  CHANGE  THE  PRECISION,  MODIFY  THE  DIMENSIONS  IN  THE 
C DECLARE  STATEMENTS  IN  THE  AUGMENT  DESCRIPTION  DECK  - 
C THE  DIMENSION  FOR  TYPE  MULTIPLE  SHOULD  BE  T+2  AND 
C FOR  TYPE  MULTIPAK  SHOULD  BE  INT  ((T+3)/2). 

C SEE  COMMENTS  IN  ROUTINE  MPSET  FOR  THE  NUMBER  OF  MP 
C DIGITS  REQUIRED  TO  GIVE  THE  EQUIVALENT  OF  ANY  DESIRED 
C NUMBER  Ot’  DECIMAL  PLACES. 

C •••  ON  SOME  SYSTEMS  A DECLARATION  OF  BLANK  COMMON  IN  THE  MAIN 
C •«  PROGRAM  MAY  BE  NECESSARY.  IF  SO,  DECLARE 
C •••  COMMON  MPWORK(301) 

C •••  OR,  MORE  GENERALLY, 

C •••  COMMON  MPWORK(MXR+5) 

C •••  IN  THE  MAIN  PROGRAM. 

RETURN 

END 


MP032825 
MP032827 
MP032829 
MP032831 
MP032833 
MP032835 
MP032837 
MP032839 
MP032841 
MP032843 
MP032845 
MP032847 
MP032849 
MP032851 
MP032853 
MP032855 
MP032857 
MP032859 
MP032861 
MP032863 
MP032865 
MP032867 
MP032869 
MP03287 1 
MP032873 
MP032875 
MPO 32877 
MPO 32879 
MPO 32881 
MP032883 
MP032885 
MP032887 
MPO 32 889 
MP032891 
MP032893 
MP032895 
MP032897 


C $$  «««  MPIO  »*•*»• 

SUBROUTINE  MPIO  (C,  N,  UNIT,  IFORM,  ERR) 

C IF  UNIT  .GT.  0 WRITES  C(1),  ...  , C(N)  IN  FORMAT  IFORM 
C IF  UNIT  .LE.  0 READS  C(1),  ...  , C(N)  IN  FORMAT  IFORM 
C IN  BOTH  CASES  USES  LOGICAL  UNIT  IABS(UNIT). 

C ERR  IS  RETURNED  AS  TRUE  IF  N NON-POSITIVE,  OTHERWISE  FALSE. 

C WE  WOULD  LIKE  TO  RETURN  ERR  AS  TRUE  IF  READ/WRITE  ERROR  DETECTED, 
C BUT  THIS  CAN  NOT  BE  DONE  WITH  ANSI  STANDARD  FORTRAN  (1966). 

C •••  UNIVAC  ASCII  FORTRAN  (FTN  5R1AE)  DOES  NOT  WORK  IF  IFORM 
C •••  IS  DECLARED  WITH  DIMENSION  1.  MOST  FORTRANS  DO  THOUGH. 
INTEGER  C(N) , UNIT,  IF0RM(20) 

LOGICAL  ERR 
ERR  = (N.LE.O) 

IF  (ERR)  RETURN 
IU  = IABS(UNIT) 

IF  (UNIT.GT.O)  WRITE  (IU,  IFORM)  C 
IF  (UNIT.LE.C)  READ  (IU,  IFORM)  C 
RETURN 


MP032921 
MP032923 
MP032925 
MPO 32927 
MP032929 
MP032931 
MP032933 
MP032935 
MP032937 
MPO 32939 
MP032941 
MP032943 
MP032945 
MP032947 
MP032949 
MP032951 
MP032953 
MP032955 


END 


MP032957 


C $$  •••••«  MPKSTR  MP032961 

SUBROUTINE  MPKSTR  (X,  Y)  MP032963 

C SETS  Y s X FOR  PACKED  MP  NUMBERS  X AND  Y.  MP032965 

C ASSUMES  SAME  PACKED  FORMAT  AS  MPPACK  AND  MPUNPK.  MP032967 

COMMON  B,  T,  M,  LUN,  MXR,  R MP032969 

INTEGER  B,  T,  R(1),  X(2),  Y(2)  MP032972 

Y(2)  = X(2)  MP032973 

C CHECK  FOR  ZERO  MP032975 

IF  (Y(2).EQ.O)  RETURN  MP032977 

C HERE  X NONZERO  SO  MOVE  PACKED  NUMBER  MP032979 

N = (T+3)/2  MP032981 

DO  10  I = 1,  N MP032983 

10  Y(I)  = X(I)  MP032985 

RETURN  MP032987 

END  MP032989 

C $$  ••••**  MPLE  •****•  MP033001 

LOGICAL  FUNCTION  MPLE  (X,  Y)  MP033003 

C RETURNS  LOGICAL  VALUE  OF  (X  .LE.  Y)  FOR  MP  X AND  Y.  MP033005 

INTEGER  X(1),  Y ( 1 ) MP033007 

MPLE  = (MPCOMP(X, Y)  .LE.  0)  MP033009 

RETURN  MP033011 

END  MP033013 

C M •*••••  MPLT  •••**•  MP037281 

LOGICAL  FUNCTION  MPLT  (X,  Y)  MP037283 

C RETURNS  LOGICAL  VALUE  OF  (X  .LT.  Y)  FOR  MP  X AND  Y.  MP037285 

INTEGER  X(1),  Y( 1 ) MP037287 

MPLT  = (MPCOMP(X,Y)  .LT.  0)  MP037289 

RETURN  MP037291 

END  MP037293 

C $$  MPMEXA  ••»•••  MP038051 

FUNCTION  MPMEXA  (X)  MP038053 

C RETURNS  THE  MAXIMUM  ALLOWABLE  EXPONENT  OF  MP  NUMBERS  (THE  THIRD  MP038055 

C WORD  OF  COMMON).  X IS  A DUMMY  MP  ARGUMENT.  MP038057 

COMMON  B,  T,  M,  LUN,  MXR,  R MP038059 

INTEGER  B,  T,  R(1),  X(1)  MPO38O6I 

MPMEXA  = M MP038063 

RETURN  MP038065 

END  MP038067 

C $$  #•••••  MPMEXB  •••*••  MP038071 

SUBROUTINE  MPMEXB  (I,  X)  MP038073 

C SETS  THE  MAXIMUM  ALLOWABLE  EXPONENT  OF  MP  NUMBERS  (I.E.  THE  MP038075 

C THIRD  WORD  OF  COJtiON)  TO  I.  MP038077 

C I SHOULD  BE  GREATER  THAN  T,  AND  4«I  SHOULD  BE  REPRESENTABLE  MP038079 

C AS  A SINGLE-PRECISION  INTEGER.  MPO38O8I 

C X IS  A DUMMY  MP  ARGUMENT  (AUGMENT  EXPECTS  ONE).  MP038083 

COMMON  B,  T,  M,  LUN,  MXR,  R MP038085 

INTEGER  B,  T,  R(1),  X(1)  MP038087 

M = I MP038089 

C CHECK  LEGALITY  OF  M.  IF  TOO  LARGE,  4»M  MAY  OVERFLOW  AND  TEST  .LE.  0 MP038091 

IF  ( (M.GT.T) .AND. ( (4*M) .GT.O) ) RETURN  MP038093 
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WRITE  (LUN,  10)  MP038095 

10  FORMAT  (44H  ATTEMPT  TO  SET  ILLEGAL  MAXIMUM  EXPONENT,  MP038097 

$ 22H  IN  CALL  TO  MPMEXB  »••)  MP038099 

CALL  MPERR  MPO38IOI 

RETURN  MP038103 

END  MP038105 

C $$  «««  MPNE  »•»«»  MP040461 

LOGICAL  FUNCTION  MPNE  (X,  Y)  MP040463 

C RETURNS  LOGICAL  VALUE  OF  (X  .NE.  Y)  FOR  MP  X AND  Y.  MP040465 

INTEGER  X(1),  Y(1)  MP040467 

MPNE  = (MPCOMP(X, Y)  .NE.  0)  MP040469 

RETURN  MP040471 

END  MP040473 

C $$  ••••••  MPOUTF  ******  MP041781 

SUBROUTINE  MPOUTF  (X,  P,  N,  IFORM,  ERR)  MP041783 

C WRITES  MP  NUMBER  X ON  LOGICAL  UNIT  LUN  (FOURTH  WORD  OF  COMMON)  MP041785 

C IN  FORMAT  IFORM  AFTER  CONVERTING  TO  FP.N  DECIMAL  REPRESENTATION  MP041787 

C USING  ROUTINE  MPOUT.  FOR  FURTHER  DETAILS  SEE  COMMENTS  IN  MPOUT.  MP041789 

C IFORM  SHOULD  CONTAIN  A FORMAT  WHICH  ALLOWS  FOR  OUTPUT  OF  P MP041791 

C WORDS  IN  A1  FORMAT,  PLUS  ANY  DESIRED  HEADINGS,  SPACING  ETC.  MP041793 

C E.G.  24H(8H1HEADING/( 1 IX, 100A1))  MP041795 

C ERR  RETURNED  AS  TRUE  IF  P NOT  POSITIVE,  OTHERWISE  FALSE.  MP041797 

C SPACE  REQUIRED  3T+P+1 1 WORDS.  MP041799 

COMMON  B,  T,  M,  LUN,  MXR,  R MP041801 

INTEGER  B,  T,  R(1),  X(1),  IFORM(I),  P MP041803 

LOGICAL  ERR  MP041805 

ERR  = .TRUE.  MP041807 

C RETURN  WITH  ERROR  FUG  SET  IF  OUTPUT  FIELD  WIDTH  P NOT  POSITIVE  MP041809 

IF  (P.LE.O)  RETURN  MP041811 

C CHECK  THAT  ENOUGH  SPACE  IS  AVAILABLE  MP041813 

CALL  MPCHK  (3,  P+11)  MP041815 

12  = 3*T  + 12  MP041817 

C CONVERT  X TO  DECIMAL  FORM  MP041819 

CALL  MPOUT  (X,  R(I2),  P,  N)  MP041821 

C AND  WRITE  ON  UNIT  LUN  WITH  FORMAT  IFORM  MP041823 

CALL  MPIO  (R(I2),  P,  LUN,  IFORM,  ERR)  MP041825 

RETURN  MP041827 

END  MP041829 

C $$  »•••«  MPSIGA  ••••••  MP048741 

FUNCTION  MPSIGA  (X)  MP048743 

C RETURNS  SIGN  OF  MP  NUMBER  X MP048745 

INTEGER  X( 1 ) MP048747 

MPSIGA  = X(1)  MP048749 

RETURN  MP048751 

END  MP048753 

C $$  ••••••  MPSIGB  ••••••  MP048761 

SUBROUTINE  MPSIGB  (I,  X)  MP048763 

C SETS  SIGN  OF  MP  NUMBER  X TO  I.  MP048765 

C I SHOULD  BE  0,  +1  OR  -1.  MP048767 

C EXPONENT  AND  DIGITS  OF  X ARE  UNCHANGED,  MP048769 

C BUT  RESULT  MUST  BE  A VALID  MP  NUMBER.  MP046771 

COMMON  B,  T,  M,  LUN,  MXR,  R MP048773 
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INTEGER  B,  T,  R(1),  X(3)  MP048775 

X(  1 ) = I MP048777 

C CHECK  FOR  VALID  SIGN  . MP048779 

IF  (IABS(I).LE.I)  GO  TO  20  MP048781 

WRITE  (LUN,  10)  MP048783 

10  FORMAT  (39H  •••  INVALID  SIGN  IN  CALL  TO  MPSIGB  •*•)  MP048785 

GO  TO  40  MP048787 

C RETURN  IF  X ZERO  MP048789 

20  IF  (I.EQ.O)  RETURN  MP048791 

C CHECK  FOR  VALID  EXPONENT  AND  LEADING  DIGIT  MP048793 

IF  ((IABS(X(2)).LE.M).AND.(X(3).GT.0).AND.(X(3).LT.B))  RETURN  MP048795 

WRITE  (LUN,  30)  MP048797 

30  FORMAT  (48H  •••  X NOT  VALID  MP  NUMBER  IN  CALL  TO  MPSIGB  •*•)  MP048799 

40  CALL  MPERR  MP048801 

X(1)  = 0 MP048803 

RETURN  MP048805 

END  MP048807 


C $$  ••••••  MPUPK  •••••• 

SUBROUTINE  MPUPK  (SOURCE,  DEST,  LDEST,  LFIELD) 

C 

C HHiHHimimimiiii 

C «•  MACHINE  DEPENDENT  ••• 

C *•••••••••**••*•••••••*•• 

C 

C MACHINE-DEPENDENT  STATEMENTS  ARE  SURROUNDED  BY  C •••  LINES 

C »*• 

C THIS  IS  UNIVAC  1100,  FORTRAN  V VERSION. 

C »•* 

C THIS  SUBROUTINE  UNPACKS  A PACKED  HOLLERITH  STRING  (SOURCE) 

C PLACING  ONE  CHARACTER  PER  WORD  IN  THE  ARRAY  DEST  (AS  IF  READ  IN 
C A 1 FORMAT).  IT  CONTINUES  UNPACKING  UNTIL  IT  FINDS  A SENTINEL  ($) 
C OR  UNTIL  IT  FINDS  A COMPILER  GENERATED  SENTINEL  (IF  SO 
C IMPLEMENTED)  OR  UNTIL  IT  HAS  FILLED  LDEST  WORDS  OF  THE 
C ARRAY  DEST.  THE  LENGTH  OF  THE  UNPACKED  STRING  IS  RETURNED 
C IN  LFIELD.  THUS  0 .LE.  LFIELD  .LE.  LDEST. 

INTEGER  SOURCE(I),  DEST(1),  BUNKS,  TEMP 
DATA  BLANKS  /1H  /,  1ST  /1H$/ 

C NK  IS  THE  NUMBER  OF  CHARACTERS  PER  WORD 
C AND  ISTC  IS  THE  COMPILER-GENERATED  SENTINEL  (IF  ANY) 

C ••• 

DATA  NK  /6/,  ISTC  /O/ 

C ••• 

TEMP  s BUNKS 
LD  = LDEST 
LFIELD  s 0 
IF  (LD.LE.O)  RETURN 
DO  10  K s 1,  LD 
I s LFIELD/NK  ♦ 1 

C GET  NEXT  WORD  (CONTAINING  NK  CHARACTERS)  AND 
C CHECK  FOR  COMPILER-GENERATED  END-OF-STRING  SENTINEL 
IF  (SOURCE(I)  .EQ.  ISTC)  RETURN 
C MOVE  (M0D(LFIELD,NK)+1 )-TH  CHARACTER  OF  SOURCE(I)  TO 
C FIRST  (I.E.  LEFTMOST)  CHARACTER  POSITION  OF  TEMP 
C »H 

FLD  (0,  6,  TEMP)  3 FLD  (6»MOD( LFIELD, 6) , 6,  SOURCE(I)) 


« 


MP052341 

MP052343 

MP052345 

MP052347 

MP052349 

MP052351 

MP052353 

MP052355 

MP052357 

MP052359 

MP052361 

MP052363 

MP052365 

MP052367 

MP052369 

MP052371 

MP052373 

MP052375 

MP052377 

MP052379 

MP052381 

MP052383 

MP052385 

MP052387 

MP052389 

MP052391 

MP052393 

MP052395 

MP052397 

MP052399 

MP052401 

MP052403 

MP052405 

MP052407 

MP052409 

MP05241 1 

MP052413 

MP052415 
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C CHECK  FOR  END-OF-STRING  SENTINEL 
IF  (TEMP  .EQ.  1ST)  RETURN 
LFIELD  = K 
10  DEST(K)  = TEMP 
RETURN 
END 

SUBROUTINE  MPUPK  (SOURCE,  DEST,  LDEST , LFIELD) 

C 

C mmiitmiHiiiHmii 

C **•  MACHINE  DEPENDENT  »• 

C HHiimHiHiimiiim 

C 

C MACHINE-DEPENDENT  STATEMENTS  ARE  SURROUNDED  BY  C LINES 

C »*• 

C THIS  IS  UNIVAC  1100,  ASCII  FORTRAN  VERSION. 

C *« 

C THIS  SUBROUTINE  UNPACKS  A PACKED  HOLLERITH  STRING  (SOURCE) 

C PLACING  ONE  CHARACTER  PER  WORD  IN  THE  ARRAY  DEST  (AS  IF  READ  IN 
C A1  FORMAT).  IT  CONTINUES  UNPACKING  UNTIL  IT  FINDS  A SENTINEL  ($) 
C OR  UNTIL  IT  FINDS  A COMPILER  GENERATED  SENTINEL  (IF  SO 
C IMPLEMENTED)  OR  UNTIL  IT  HAS  FILLED  LDEST  WORDS  OF  THE 
C ARRAY  DEST.  THE  LENGTH  OF  THE  UNPACKED  STRING  IS  RETURNED 
C IN  LFIELD.  THUS  0 .LE.  LFIELD  .LE.  LDEST. 

INTEGER  SOURCE(I),  DEST(1),  BUNKS,  TEMP 
DATA  BLANKS  /1H  /,  1ST  /1H$/ 

C NK  IS  THE  NUMBER  OF  CHARACTERS  PER  WORD 
C AND  ISTC  IS  THE  COMPILER-GENERATED  SENTINEL  (IF  ANY) 

C ••• 

DATA  NK  /4/ , ISTC  /0/ 

C «•* 

TEMP  = BUNKS 
LD  = LDEST 
LFIELD  = 0 
IF  (LD.LE.O)  RETURN 
DO  10  K r 1,  LD 
I = LFIELD/NK  ♦ 1 

C GET  NEXT  WORD  (CONTAINING  NK  CHARACTERS)  AND 
C CHECK  FOR  COMPILER-GENERATED  END-OF-STRING  SENTINEL 
IF  (SOURCE(I)  .EQ.  ISTC)  RETURN 
C MOVE  (MOD( LFIELD, NK)+1 )-TH  CHARACTER  OF  SOURCE(I)  TO 
C FIRST  (I.E.  LEFTMOST)  CHARACTER  POSITION  OF  TEMP 
C *« 

||§  BITS  (TEMP,  1,  9)  = BITS  (SOURCE(I),  9»MOD(LFIELD,4)+1 , 9) 

C CHECK  FOR  END-OF-STRING  SENTINEL 
IF  (TEMP  .EQ.  1ST)  RETURN 
LFIELD  = K 
10  DEST(K)  = TEMP 
RETURN 
END 

C $$  MPUPK  •••••• 

SUBROUTINE  MPUPK  (SOURCE,  DEST,  LDEST,  LFIELD) 


MP052417 

MP052419 

MP052421 

MP052423 

MP052425 

MP052427 

MP052429 


MP052341 

MP052343 

MP052345 
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C IIIHIHHHHIflHHIHI  MP052347 

C •«  MACHINE  DEPENDENT  •••  MP052349 

C IHIHIIIHIHHIHtllHI  MP052351 

C MP052353 

C MACHINE-DEPENDENT  STATEMENTS  ARE  SURROUNDED  BY  C LINES  MP052355 

C MP052357 

C THIS  IS  IBM  360  FORTRAN  G OR  H VERSION 

C »«  MP052361 

C THIS  SUBROUTINE  UNPACKS  A PACKED  HOLLERITH  STRING  (SOURCE)  MP052363 

C PLACING  ONE  CHARACTER  PER  WORD  IN  THE  ARRAY  DEST  (AS  IF  READ  IN  MP052365 

C A1  FORMAT).  IT  CONTINUES  UNPACKING  UNTIL  IT  FINDS  A SENTINEL  ($)  MP052367 

C OR  UNTIL  IT  FINDS  A COMPILER  GENERATED  SENTINEL  (IF  SO  MP052369 

C IMPLEMENTED)  OR  UNTIL  IT  HAS  FILLED  LDEST  WORDS  OF  THE  MP052371 

C ARRAY  DEST.  THE  LENGTH  OF  THE  UNPACKED  STRING  IS  RETURNED  MP052373 

C IN  LFIELD.  THUS  0 .LE.  LFIELD  .LE.  LDEST.  MP052375 


INTEGER  DEST ( 1 ) , BLANKS,  TEMP 
C 

LOGICAL* 1 SOURCE ( 1 ) , TC(4) 
EQUIVALENCE  (TC,  TEMP) 

C ••• 


DATA  BLANKS  /1H  /,  1ST  /1H$/  MP052379 

C NK  IS  THE  NUMBER  OF  CHARACTERS  PER  WORD  MP052381 

C AND  ISTC  IS  THE  COMPILER-GENERATED  SENTINEL  (IF  ANY)  MP052383 

C •••  MP052385 

DATA  NK  /4/,  ISTC  /O/ 

C •••  MP052389 

TEMP  = BUNKS  MP052391 

LD  = LDEST  MP052393 

LFIELD  = 0 MP052395 

IF  (LD.LE.O)  RETURN  MP052397 

DO  10  K = 1,  LD  MP052399 

C •••  MP052413 

TC( 1 ) = SOURCE  (K) 

C •••  MP052417 

C CHECK  FOR  END-OF-STRING  SENTINEL  MP052419 

IF  (TEMP  .EQ.  1ST)  RETURN  MP052421 

LFIELD  = K MP052423 

10  DEST(K)  = TEMP  MP052425 

RETURN  MP052427 

END  MP052429 
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